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Introduction 


Data  obtained  from  high  altitude  satellites  have  shown  that  the  atmo- 
spheric density  has  a strong  diurnal  variation,  maximizing  near  the  subsolar 
point  at  around  14  LST.  For  the  most  part,  atmospheric  models  which  have 
been  based  principally  on  data  gathered  above  about  250  km,  where  the  diurnal 
character  is  very  pronounced,  have  extrapolated  this  behavior  to  lower  alti- 
tudes though  with  decreasing  amplitude.  Some  indication  that  such  a straight- 
forward extrapolation  may  not  be  correct  has  been  given  by  the  San  Marco  3 
satellite  data  presented  by  Newton  et  al.  (1975).  These  investigators  found  at 
220  km  two  distinct  peaks  of  equal  amplitude  at  1 1 LST  and  18.5  LST.  Newton 
et  al.  attributed  this  behavior  to  the  different  phase  in  the  diurnal  response  of 
individual  constituents  O and  N^*  That  the  number  density  of  each  atmospheric 
constituent  peaks  at  a different  time  of  day  has  also  been  deducted  by  Mayr 
et  al.  (1974)  from  analysis  of  OGO-6  mass  spectrometer  data. 

We  have  used  density  data  from  a cold  cathode  gauge  on  board  the  equa- 
torial satellite  Atmosphere  Explorer  E to  investigate  the  diurnal  behavior  of 
the  atmosphere  from  300  km  down  to  140  km.  The  data  are  used  to  determine 
the  phase  and  relative  amplitude  of  the  diurnal  and  semidiurnal  components  of 
the  local  time  variation  as  a function  of  altitude.  These  parameters  are  valuable 
as  an  indication  of  the  relative  magnitude  of  energy  sources  for  the  lower  thermo- 
sphere since  a diurnal  component  may  be  expected  to  result  from  local  heating 


source  would  also  depend,  of  course,  on  such  detcdls  as  phase  and  phase 
progression  and  the  altitude -dependence  of  the  wave  amplitude. 
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Discussion 


The  data  used  in  this  study  consist  of  5062  density  measurements  obtained 
using  the  cold  cathode  ion  gauge  known  as  Pressure  Sensor  A (PSA)  on  board 
the  Atmosphere  Explorer  E (AE-E)  satellite.  The  instrument  has  been  de- 
scribed by  Rice  et  al.  (19  73),  and  the  data  reduction  procedure  has  been  pre- 
sented by  Rice  and  Sharp  (1977).  AE-E  is  in  a near -equatorial  orbit,  with  an 
inclination  of  19.6°.  The  data  used  in  this  study  were  gathered  between 
December  1975  and  October  1976  when  AE-E  was  in  an  eccentric  orbit  (140  km 
X 275  7 km).  The  data  set  was  further  limited  by  including  only  data  obtained 
during  periods  of  oriented  (non -spinning)  operation.  Density  and  ephemeris 
data  at  fixed  altitude  levels  at  5 km  intervals  (140-300  km)  were  tabulated  for 
use  in  the  present  study.  In  order  to  minimize  the  complicating  and  incom- 
pletely modeled  effects  of  geomagnetic  activity,  the  data  base  excludes  any 
data  for  which  Kp  > 5 at  any  time  during  the  period  six  hours  prior  to  the  time 
of  measurement. 

In  order  to  isolate  the  loccil  time  variation  in  the  lower  thermosphere, 
all  the  other  known  contributions  to  the  variability  of  the  density  were  sup- 
pressed using  the  Jacchia  (1971)  model,  evaluated  using  all  appropriate  input 
parameters,  but  with  the  local  time  set  arbitrarily  at  midnight.  The  ratio  R 
of  the  measured  density  to  this  model  value  was  then  averaged  over  each  two- 
hour  interval  of  local  time.  The  results  of  this  analysis  are  shown  in  Figure  1. 
Each  curve  represents  the  average  of  this  ratio  over  the  five  altitude  levels 
in  a 20  km  band  centered  at  the  designated  altitude.  (For  example,  the  290  km 
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curve  includes  data  from  the  280-300  km  range.  ) The  curves  thus  represent 
the  qualitative  altitude  dependence  of  the  local  time  variation  which  occurs 
near  the  equator.  There  is  a systematic  transition  from  the  clearly  diurnal 
behavior  at  290  km  with  a well-defined  maximum  near  14  hours  local  time 
and  a broader  minimum  near  midnight  to  the  semidiurnal  behavior  at  140  km. 

In  the  transition  region  between  these  altitudes  the  daily  variation  is  some- 
what less  well-defined  due  to  the  fact  that  neither  the  diurnal  nor  the  semi- 
diurnal mode  is  clearly  dominant. 

In  order  to  study  the  altitude  dependence  of  the  daily  density  variation  in 
greater  detail,  we  have  performed  a two -component  harmonic  analysis  of  the 
ratio  R after  dividing  by  the  mean  value.  The  amplitude  and  phase  of  the  di- 
urnal and  semidiurnal  terms  were  calculated  for  each  5 km  altitude  step  in 
the  140-300  km  data  base.  Initially  a three -component  fit  was  made,  but  the 
terdiurnal  component  showed  no  phase  continuity  with  altitude,  and  the  ampli- 
tude at  all  points  was  too  small  to  be  statistically  significant.  Furthermore, 
neglect  of  the  third  harmonic  in  later  analyses  did  not  significantly  affect  either 
the  amplitude  or  phase  of  the  resultant  diurnal  and  semidiurnal  terms.  A 
somewhat  similar  analysis  using  data  from  the  San  Marco  3 satellite  at  higher 
altitudes  has  been  described  by  Broglio  et  al.  (19  76). 

In  Figures  2 and  3 we  present  the  amplitude  and  phase,  respectively, 
determined  in  each  of  these  decompositions.  The  amplitudes  and  phases  of 
the  two  modes  are  shown  as  functions  of  altitude,  along  with  the  correspondi.'^g 
predictions  of  two  widely  used  empirical  models  of  the  thermosphere,  the 
Jacchia  (1971)  model  and  the  Mass  Spectrometer -Incoherent  Scatter  (MSIS) 
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of  fit  at  each  altitude;  for  clarity  they  are  only  given  every  other  altitude 
increment.  Also  plotted  are  the  results  of  performing  a similar  decompositic 
on  the  density  values  calculated  by  means  of  the  Jacchia  (1971)  and  MSIS 
thermospheric  models. 


al  and  sennidiurnal  components  of  the  local -time 
nding  to  the  amplitudes  in  Fig.  2. 


i 

I 

model  of  Hedin  et  al.  (1977).  In  calculating  the  amplitude  and  phase  predic-  i 

tions  of  the  empirical  models,  model  densities  are  determined  for  every  point  ! 

in  our  AE-E  data  set  and  the  resultant  values  are  then  handled  in  the  same  ! 

I 

manner  as  with  the  actual  measured  data. 

1 

Comparison  of  the  amplitudes  derived  from  the  data  with  those  of  the 
empirical  models  leads  to  several  conclusions.  For  both  the  data  and  the 

models,  the  amplitude  of  the  diurnal  mode  decreases  rapidly  with  decreasing  i 

altitude.  However,  at  any  given  altitude  the  amplitude  of  the  diurnal  model  i 

i 

predicted  by  the  Jacchia  (19  71)  model  exceeds  that  actually  measured  by  more  ^ 

than  a factor  of  two.  The  altitude  dependence  of  the  amplitude  of  the  MSIS  j 

1 

diurnal  mode  closely  resembles  that  of  the  AE-E  data.  Although  the  ampli- 

• I 

tude  is  uniformly  less  than  the  measured  values,  it  is  only  slightly  beyond  : 

the  uncertainty  limits  for  the  fits.  The  amplitudes  of  the  semidiurnal  mode  j 

calculated  from  the  data  and  both  models  are  comparable  at  the  higher  alti- 
tudes treated  here  (>240  km).  However,  this  overall  agreement  breaks  down 
below  ~200  km,  where  the  data  and  the  MSIS  model  indicate  an  increase  with 
decreasing  altitude,  whereas  the  semidiurnal  mode  in  the  Jacchia  (19  71)  model 
continues  to  decrease  monotonic  ally.  The  measured  semidiurnal  mode  sur- 
passes the  diurnal  mode  in  amplitude  at  ~ 180  km  and  is  dominant  at  all  lower 
altitudes.  The  essentially  identical  behavior  of  the  semidiurnal  mode  in  the 
MSIS  model  and  in  our  AE-E  data  base  extends  to  such  details  as  the  altitude 
at  which  minimum  amplitude  is  seen  (~200  km).  The  agreement  of  model  and 
data  is  even  more  striking  in  the  phase  plot.  Figure  3.  The  phase  of  the  semi- 
diurnal mode,  which  at  high  altitudes  is  almost  constant  between  2 and  3 hours 
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LST,  makes  near  200  km  a nearly  discontinuous  transition  to  a phase  angle 
between  8 and  10  hours  LST.  The  Jacchia  (1971)  model  extrapolates  the  be- 
havior at  higher  altitudes  downward  to  the  lower  thermosphere  as  well. 

The  measured  phase  of  the  diurnal  mode  shows  good  agreement  with 
the  models  at  high  altitudes,  but  considerable  disagreement  at  low  altitudes. 
The  phase  in  the  measured  data  progresses  from  near  13  hours  LST  at  300  km 
toward  earlier  times,  reaching  6 hours  LST  at  140  km.  The  phase  of  the 
MSIS  model,  on  the  other  hand,  though  starting  at  a similar  point  at  high  alti- 
tude, progresses  toward  later  times,  reaching  2.5  hours  LST  at  140  km.  The 
Jacchia  (1971)  model  shows  a diurnal  phase  which  remains  constant  at  14  hours 
LST.  This  divergence  can  probably  be  attributed  to  the  manner  in  which  the 
data  used  in  constructing  the  models  was  obtained.  The  data  base  for  the 
Jacchia  (1971)  model  conteiined  no  data  below  an  altitude  of  ~260  km;  the  con- 
stant phase  results  from  extrapolating  downwards  the  behavior  of  the  upper 
thermosphere  and,  hence,  good  agreement  in  the  lower  thermosphere  might 
not  be  expected.  On  the  other  hand,  the  low  altitude  data  input  to  the  MSIS 
model  is  almost  exclusively  from  the  Atmosphere  Explorer  C satellite.  We 
believe  that  the  behavior  of  the  phase  of  the  diurnal  mode  in  the  MSIS  model 
may  be  an  artifact  arising  from  the  strong  coupling  in  the  AE-C  data  base  of 
local  time  with  other  parameters  such  as  latitude  and  season.  Because  the 
amplitude  of  the  diurnal  mode  is  small  at  low  thermospheric  altitudes,  such 
behavior  has  little  effect  on  the  density  there.  To  test  this  hypothesis  we 
performed  the  above  analysis  using  our  own  ion  gauge  data  from  the  AE-C 
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satellite  and  obtained  results  which  were  qualitatively  sinailar  to  the  MSIS 
predictions  with  respect  to  the  diurnal  phase  progression  at  low  altitude. 

The  MSIS  model  semidiurnal  mode  agrees  quite  well  in  phase  with  the 
data  showing  the  sharp  phase  transition  associated  with  the  change  from  EUV 
to  tidal  domination  below  ~Z00  km.  The  Jacchia  (1971)  model  fails  to  repro- 
duce this  behavior  and  disagrees  substantially  in  phase  at  the  lower  altitudes 
studied  here.  This  suggests  that  the  MSIS  model  is  more  reliable  for  studies 
involving  lower  thermospheric  chemistry  and  dynamics  than  is  the  Jacchia 
(1971)  model.  In  this  connection  it  should  be  noted  that  the  MSIS  model  ex- 
hibits the  semidiurnal  mode  as  dominant  only  at  low  latitude;  at  140  km  the 
diurnal  mode  is  dominant  above  about  30°  latitude.  Our  data  from  AE-E 
cannot  confirm  or  deny  this  behavior,  however. 

The  general  behavior  of  a dominant  diurnal  mode  at  high  altitudes  and 
semidiurnal  mode  at  low  altitudes  is  consistent  with  the  following  picture. 

At  high  altitudes  the  characteristics  of  thermosphere  are  determined  largely 
by  the  day -night  variations  in  EUV  heating,  leading  to  a diurnal  variation  whose 
amplitude  decreases  with  decreasing  altitude  because  the  heating  rate  per  unit 
mass  decreases  with  decreasing  altitude.  The  behavior  at  low  thermosphere 
altitudes  arises  from  the  propagation  of  tidal  waves  from  lower  regions  of  the 
atmosphere.  These  tides  are  known  to  be  predominantly  semidiurnal  in  char- 
acter (Lindzen  et  al.  1976;  Wcind,  1976)  and  their  amplitudes  decrease  with 
increasing  altitude  in  the  thermosphere. 

Although  comparison  of  phase  in  density  data  with  that  in  temperature 
must  be  approached  with  a good  deal  of  caution,  it  is  worth  pointing  out  the 
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qualitative  similarity  of  the  altitude  dependence  of  the  density  phase 
measurements  reported  here  and  the  temperature  phase  calculated  for  the 
solar  semidiurnal  tide  by  Lindzen  and  Hong  (1974)  for  equinoctial  conditions 
near  the  equator,  including  the  nearly  180°  phase  change  seen  near  200  km 
and  the  more  gradual  shift  toward  later  local  time  evident  below  160  km.  The 
calculations  of  temperature  phase  by  Lindzen  and  Hong  (their  Figure  20)  show 
this  behavior  uniformly  40-50  km  below  the  ailtitudes  at  which  it  is  seen  in  our 
density  measurements.  On  the  other  hand,  the  calculations  for  solstitial  con- 
ditions show  much  more  extreme  phase  changes  with  altitude  (720°  over  the 
same  altitude  band,  125-175  km),  and  since  our  data  base  encompasses  sol- 
stitial and  equinoctial  conditions  approximately  equally,  the  previous  noted 
qualitative  agreement  may  be  only  fortuitous.  Given  the  general  character 
of  the  altitude  dependence  of  the  semidiurnal  mode,  it  seems  to  us  reasonable 
to  ascribe  this  behavior  to  the  solar  semidiurnal  tide.  With  this  interpretation, 
the  transition  between  the  part  of  the  thermosphere  dominated  by  tides  and  the 
EUV -dominated  upper  thermosphere  occurs  near  the  equator  in  the  200-220  km 
altitude  range. 
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Conclusion 


The  Atmosphere  Explorer  E density  data  suggest  that  the  upward 
propagating  mesospheric  tidal  wave  dominates  the  local  time  behavior  near 
the  equator  in  the  lower  thermosphere  up  to  about  180  km  and  is  important  to 
at  least  200  km.  Thus  the  common  assumption  of  a predominantly  diurnal 
density  variation  below  ~ 200  km  may  lead  to  serious  errors  in  density  pre- 
dictions and  cause  errors  in  any  calculation  attempting  to  determine  the  dynam- 
ical behavior  of  the  atmosphere  below  200  km.  Of  those  atmospheric  models 
known  to  us,  only  the  MSIS  model  of  Hedin  et  al.  (1977)  reproduces  the  mea- 
sured data. 
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THE  IVAN  A,  GETTING  l^BORATORIES 


The  Laboratory  Operations  of  The  Aerospace  Corporation  is  conducting 
experimental  and  theoretical  investigations  necessary  for  the  evaluation  and 
application  of  scientific  advances  to  new  military  concepts  and  systems.  Ver> 
satility  and  flexibility  have  been  developed  to  a high  degree  by  the  laboratory 
personnel  in  dealing  with  the  many  problems  encountered  in  the  nation's  rapidly 
developing  space  and  missile  systems.  Expertise  in  the  latest  scientific  devel- 
opments is  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 

Aerophysics  Laboratory;  Launch  and  reentry  aerodynamics,  heat  trans- 
fer, reentry  physics,  chemical  kinetics,  structural  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high-power  gas  lasers. 

Chemistry  and  Physics  Laboratory;  Atmospheric  reactions  and  atmos- 
pheric  optics,  chemical  reactions  in  polluted  atmospheres,  chemical  reactions 
of  excited  species  in  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo- 
sensitive materials  and  sensors,  high  precision  laser  ranging,  and  the  appli- 
cation of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 

Electronics  Research  Laboratory;  Electromagnetic  theory,  devices,  and 
propagation  phenomena,  including  plasma  electromagnetics;  quantum  electronic j, 
lasers,  and  electro-optics;  communication  sciences,  applied  electronics,  semi- 
conducting, superconducting,  and  crystal  device  physics,  optical  and  acoustical 
imaging;  atmospheric  pollution;  millimeter  wave  and  far-infrared  technology. 

Materials  Sciences  Laboratory:  Development  of  new  materials;  metal 
matrix  cornposites  and  new  forms  of  carbon;  test  and  evaluation  of  graphite 
and  ceramics  in  reentry;  spacecraft  materials  and  electronic  components  in 
nuclear  weapons  environment:  application  of  fracture  mechanics  to  stress  cor- 
rosion and  fatigue-induced  fractures  in  structural  metals. 

Space  Sciences  Laboratory;  Atmospheric  and  ionospheric  physics,  radia- 
tion from  the  atmosphere,  density  and  composition  of  the  atmosphere,  aurorae 
and  airglow;  magnetospheric  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  waves  in  the  magnetosphere:  solar  physics,  studies  of  solar  magnetic 
fields;  space  astronomy,  x-ray  astronomy;  the  effects  of  nuclear  explosions, 
magnetic  storms,  and  solar  activity  on  the  earth's  atmosphere,  ionosphere,  and 
magnetosphere;  the  effects  of  optical,  electromagnetic,  and  particulate  radia- 
tions in  space  on  space  systems. 
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